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SUMMARY 


An  approximate  closed  form  analysis  is  used  to  determine  the 

practical  ideal  hovering  performance  of  a  circulation  control  rotor. 

The  rotor  employs  a  cambered  elliptical  airfoil  with  a  high  section 

equivalent  lift-drag  ratio.  This  feature  is  shown  to  be  critical 

for  good  hover  performance.  The  calculated  performance  is  compared 

with  results  from  a  more  exact  prediction  theory  and  shown  to  yield 

good  agreement.  Finally,  a  comparison  is  made  with  a  rotor  using  a 

standard  NACA  0012  section.  It  is  demonstrated  that  a  circulation 

* 

control  rotor  can  achieve  comparable  hover  efficiency  to  a  conven¬ 
tional  rotor  but  over  a  much  greater  thrust  range. 
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INTRODUCTION 


There  is  a  present  need  to  quantify  the  hover  efficiency  for  a 
Circulation  Control  Rotor  (CCR) .  A  methodology  will  be  established  to 
determine  the  maximum  obtainable  figure  of  merit  of  a  CCR  and  to  make  a 
relative  comparison  between  a  conventional  helicopter  rotor  system 
utilizing  a  NACA  0012  airfoil  and  a  blown  elliptical  airfoil  section. 

The  maximum  figure  of  merit  defined  in  this  report  is  in  fact  a  "prac¬ 
tical  ideal"  figure  of  merit  which  a  rotor  could  obtain  if  designed 
purely  for  hover.  The  significance  of  such  an  analysis  is  to  provide  a 
datum  by  which  more  practical  rotor  designs  may  be  compared.  The  latter 
will  generally  vary  from  the  ideal  hover  design  due  to  the  inevitable 
compromises  for  forward  flight,  autorotation,  structure  and  dynamics. 

In  order  to  quantify  the  hover  performance  closed  form  analytical 
expressions  were  developed  which  may  be  used  to  gain  a  better  under¬ 
standing  of  the  factors  which  determine  hover  efficiency.  Simple  inte¬ 
gral  expressions  are  developed  herein;  the  only  major  assumptions  being 
a  linear  chordwise  taper  ratio  and  ideal  twist. 

The  analysis  utilizes  experimental,  two  dimensional  airfoil  data 
and  associated  equivalent  lift-drag  ratios.  This  parameter  is  defined 
in  such  a  manner  as  to  include  all  rotor  power  losses  except  for  the 
induced  and  coriolis  (pumping)  power  which  may  be  computed  separately 
as  three  dimensional  contributions.  The  equivalent  lift-drag  ratio 
(Reference  1)  is  given  by 


i/»e  =  <ycd 


where  the  equivalent  drag  coefficient  is  defined  as: 

Vj  V« 

C.  —  C,  +  C  +  C  rr— 

de  dp  M.  2V.  Vj 


(1) 

(2) 


The  first  coefficient,  C.  ,  is  the  net  profile  drag  on  the  airfoil 

P 

section  (including  any  thrusting  effect  of  the  jet  sheet).  The  second 


term,  C 

U> 


is  the  gas  horsepower  required  to  develop  lift  at  the 

V 


section  expressed  as  an  equivalent  drag. 


The  last  term, 


is 


called  the  "ram  drag"  and  is  a  small,  somewhat  pessimistic  quantity 
normally  included  in  the  definition  of  equivalent  lift  to  drag  ratio  in 
order  to  permit  a  valid  comparison  with  other  airfoils  (Reference  1). 
This  term  conservatively  accounts  for  changes  of  fluid  momentum  in  a  real 


system  amounting  to 


1 

r 


—  dx,  where  is  the  location  section 


o 

velocity.  In  the  present  analysis  of  a  rotor  blade  it  may  be  interpreted 
as  a  loss  of  the  duct  air  momentum  attributed  to  the  intake  and  com¬ 
pressor  losses  of  total  pressure.  An  additional,  loss  term  will  be  added 
under  Coriolis  Power  (p.  10)  to  account  for  the  energy  expended  to  increase 
the  angular  momentum  of  the  air  as  it  passes  out  the  duct. 


SECTION  DATA  USED 

The  CC  airfoil  section  selected  for  the  investigation .has  extremely 

good  £/d  characteristics  for  a  blown  airfoil.  Its  detailed  characteristics 
e 

and  performance  are  described  in  Reference  1  while  a  theoretical  explanation 
of  the  section  performance  is  given  in  Reference  2. 

The  CC  section  geometry  used  in  this  analysis  is  shown  in  Figure  1 
together  with  the  NACA  0012,  a  standard  helicopter  reference  airfoil. 

The  corresponding  variations  of  l/ dg  with  lift  coefficient  are  shown  in 
Figures  2  and  3.  Also  shown  in  the  figures  are  the  assumed  variation  of 
the  envelope  of  maximum  Vde*  These  analytical  expressions  are  given  by: 

jfc/d  =  C.  C  for  C  <  C  at  j t/d  (3) 

J  *  *  max 

Vde  =  ci  +  C2/Cjfc  for  Ci  >  Ci  at  A/de  W 
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In  general,  any  airfoil  section  may  be  used  in  the  analysis  which 
can  be  approximated  by  the  above  expression.  The  numerical  values  of 
these  coefficients  for  the  two  airfoils  considered  here  are: 


Circulation  Control  Ellipse 


c2  =  96 
c3  -  90 

Cj  =  1.0  at  j«/de 


NACA  0012 
C  =  -12 

C2  =  1*3.5 

c3  =  61.3 

C  =  0.75  at  x/d 

Cmax 


The  NACA  0012  airfoil  data  are  from  Reference  3  for  high  Reynolds 
number  and  standard  roughness.  These  data  are  believed  to  be  a  fairly 
accurate  representation  of  an  actual  rotor  section  and  were  found  to 
yield  very  reasonable  values  of  figure  of  merit,  comparable  to  actual 
test  values  (see:  CORRELATION  WITH  EXPERIMENTAL  DATA).  The  NACA  0012 
does  not  represent  an  optimum  conventional  section  but  Is  used  only  as 
a  convenient  datum  line.  It  should  also  be  noted  that  the  airfoil 
characteristics  are  necessarily  approximate  because  at  best  they  only 
represent  the  mean  effect  of  such  variables  as  Mach  number  (assumed 
incompressible),  Reynolds  number,  radial  flow,  turbulence,  and  manufac¬ 
turing  tolerances.  The  A/de  number  approximations  chosen  were  selected 
with  these  considerations  in  mind. 

Clearly  the  CC  airfoil  exhibits  a  much  greater  efficiency  than  the 
reference  airfoil.  To  a  large  degree  this  is  due  to  turbulent  boundary 
layer  growth  on  the  conventional  airfoil  which  contributes  to  both  the 
friction  drag  and  form  drag  and  virtually  halves  the  smooth  airfoil 
performance.  Although  the  chordal  Reynolds  number  is  usually  well  within 
the  range  for  laminar  flow  on  a  rotor  blade  in  hover,  boundary  layer 
transition  normally  occurs  near  the  leading  edge  of  conventional  airfoils 
with  the  notable  exception  of  small  blade  pitch  angles.  This  transition 
is  due  to  the  strong  adverse  pressure  gradients  associated  with  operation 
at  positive  angles  of  attack.  On  the  other  hand,  a  CC  airfoil  operates 
most  efficiently  n§ar  zero  angle  of  attack  with  very  favorable  pressure 
gradients  up  to  lift  coefficient  of  nominally  1.5.  Therefore  partial 
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laminar  flow  with  a  practical  CC  rotor  construction  is  a  very  feasible 
proposition  particularly  in  view  of  a  30-50  percent  reduction  in  Reynolds 
number  by  use  of  smaller  chord  and  tip  speed. 

A  recent  series  of  airfoil  tests  conducted  since  the  completion  of 
the  present  work  have  further  validated  the  airfoil  data.  These  tests 
covered  a  larger  range  of  Reynolds  numbers  with  many  more  test  points  and 
achieved  values  of  jt/dg  as  high  as  120.  Some  of  these  results  are  also 
shown  in  Figure  2.  Based  on  this  new  data  the  CC  approximation  would 
appear  to  be  overly  conservative  even  if  the  boundary  layer  flow  were 
assumed  fully  turbulent  (i.e.,  doubling  the  skin  friction  drag  derates 
the  i/dg  by  only  ten  percent) . 

It  is  also  important  to  note  that  recent  advances  in  conventional 
airfoils  have  tended  toward  increased  use  of  camber  to  improve  the  air¬ 
foil  efficiency,  but  at  the  expense  of  increase  pitching  moments.  In 
the  CC  airfoil  these  features  are  inherent  in  the  design  and  because  the 
elastic  axis  is  near  mid  chord  they  can  be  realized  with  relatively  small 
moments.  Reference  2  discusses  some  of  these  considerations  in  greater 
detail . 

With  regard  to  the  relative  assumptions  of  the  CC  and  NACA  0012 
airfoils  it  is  logical  to  assume  that  the  NACA  representation  is  the  more 
accurate  simply  because  it  has  worked  before  in  correlating  rotor 
performance.  The  CC  characteristics  assume  approximately  a  25  percent 
reduction  in  peak  performance  opposed  to  almost  50  percent  for  the  NACA 
0012.  However,  it  is  not  logical  in  any  way  to  derate  the  CC  by  the 
same  amount  because  the  lift  and  drag  mechanisms  are  quite  different  as 
discussed  briefly  above.  It  is,  in  any  event,  very  gratifying  that  the 
agreement  between  experiment  and  theory  for  the  NACA  0012  and  between 
the  present  theory  and  detailed  theory  for  the  CCR  are  both  in  good 
agreement . 


ANALYTICAL  APPROACH 

The  objective  is  to  define  an  optimum  envelope  of  figure  of  merit 
versus  thrust  coefficient.  This  envelope  does  not  represent  any  single 
design  but  rather  an  upper  limit  for  all  possible  configurations.  It 
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represents  that  maximum  which  could  ideally  be  achieved  by  a  given  design 
at  one  value  of  thrust  coefficient.  To  generate  this  envelope,  the  blade 
twist  is  varied  so  that  at  each  thrust  coefficient  the  blade  section 
works  at  the  angle  of  attack  for  best  efficiency  at  the  required  section 
lift  coefficient.  The  required  lift  coefficient  is  determined  to  give 
the  ideal  spanwise  loading  for  minimum  induced  power.  The  minimum 
induced  power  approach  will  yield  very  close  to  the  minimum  hover  power 
because  at  practical  rotor  thrust  coefficients  the  induced  power  amounts 
to  about  two  thirds  of  the  total  power  required.  In  theory,  under  cer¬ 
tain  specialized  conditions,  it  is  possible  to  achieve  a  somewhat  lower 
power  consumption  by  working  the  airfoil  at  a  lift  coefficient  which 
yields  smaller  profile  power  and  slightly  higher  induced  power.  However, 
such  gains  are  very  margined  and  certainly  not  within  the  accuracy  of  the 
present  work. 

In  addition,  the  effect  of  rotor  solidity  and  blade  taper  ratio  are 
evaluated.  Here  it  is  assumed  that  the  chord  varies  linearly  from  root 
to  tip.  The  combinations  of  rotor  solidity  emd  blade  taper  ratio  are 
identified  which  will  achieve  the  best  efficiency  at  a  specified  thrust 
coefficient. 

General  expressions  are  defined  for  each  component  of  rotor  power 
from  which  simplifying  assumptions  are  made  to  perform  the  actual  calcu¬ 
lations.  However,  it  will  be  shown  in  the  following  development  that 
these  assumptions  are  sufficiently  accurate  for  the  present  purpose  and 
in  fact,  agree  surprisingly  well  with  a  more  exact  analysis.  This  latter 
analysis  includes  such  effects  as:  (l)  internal  ducting  losses  as  a 
function  of  internal  duct  geometry  and  obstructions,  slot  geometry,  duct 
shape,  friction  coefficient,  duct  Mach  number,  centrifugal  compression, 
temperature  and  heat  transfer  effects;  (2)  non-uniform  induced  velocity 
distribution  computed  by  the  Lock-Goldstein  method  modified  to  account 
for  the  concentrated  tip  vortex;  (3)  variable  section  characteristics  as 
a  function  of  thickness-chord  ratio,  camber,  angle  of  attack,  blowing 
coefficient,  slot  height  to  chord  ratio,  slot  height  to  trailing  edge 
radius  ratio,  and  local  Mach  number.  Obviously  such  a  degree  of  refine¬ 
ment  is  necessary  for  a  detailed  design  study,  but  at  the  same  time,  by 
its  very  nature  it  is  quite  cumbersome  and  would  require 
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enormous  computer  time  to  make  the  present  type  of  study.  Furthermore, 
the  closed  form  expressions  developed  herein  permit  a  basic  understanding 
of  the  hover  problem,  clearly  showing  the  contribution  of  each  term.  In 
particular,  the  importance  of  the  section  equivalent  lift-drag  ratio  in 
determining  overall  rotor  efficiency  is  shown. 


HOVER  EFFICIENCY  AND  GOVERNING  EQUATIONS 

FIGURE  OF  MERIT 

The  efficiency  of  a  rotor  in  hover  is  defined  by  the  Figure  of  Merit: 


FM  = 


(5/ 


Ideal  Induced  Power 
Total  Rotor  Power 

c  3/2 

=  72c~ 


Efficient  hovering  rotors  typically  operate  at  values  of  FT4  -  0.7. 
If  the  disc  loading  and  altitude  are  known  the  power  loading,  T/HP 
(lbs/horsepower),  may  be  determined  from: 


T/HP  -  55°  FM 

T/Hp  '  yT/(’Jpg) 


(6) 


This  parameter  is  plotted  in  Figure  4  as  a  function  of  disc  loading 
and  Figure  of  Merit  for  the  standard  and  "hot  day"  conditions. 

GENERAL  POWER  EQUATIONS 

The  general  equation  for  the  total  rotor  power  required  in  flight  is: 


P=P+P+P.  +  P  +P  +P 
T  p  c  i  cor  mom  r 

where  each  of  the  power  components  are  given  below: 

R 

Profile:  P?  =  N  J V(r)  dp(r)  cos  0i(r)  dr 


(7) 


(8) 
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Compressor;  including  duct  efficiency: 


vw 


dr 


Induced: 


R 

Pi  =  N  f V(r)i(r)  sin  0^(r)  dr 


(9) 


(10) 


Coriolis: 


P 


cor 


r  m(r)  dr 


(11) 


Inlet  Momentum: 


mom 


Vi 


(12) 


Rotor  Drag: 


=  D  V 

r  oo 


(13) 


In  hover  the  last  two  terms,  Equations  (12)  and  (13)  are  zero. 

It  is  assumed  that  the  rotor  will  be  shaft  driven  so  that  the  torque  is, 
in  general,  positive.  Each  of  the  remaining  equations  (8  through  11) 
will  be  developed  into  a  closed  form  in  the  next  section. 

Induced 

The  calculation  of  induced  power  is  greatly  simplified  in  the 
present  analysis  by  the  assumption  that  the  rotor  can  generate  constant 
downwash  velocity  and  hence  develop  minimum  induced  power.  In  practice 

this  may  be  very  nearly  accomplished  by  adjusting  the  spanwise  distri¬ 
bution  of  slot  height  (Reference  4).  From  annulus  momentum  theory  the 
induced  power  is  given  by: 

R 

Pi  “  N J Vi  i£  =  J  Vi3  rdr  UM 

o 

and  for  constant  VA  an  ideal  power  coefficient  can  be  derived  as: 

C  =  0.707  C„3/2  (15) 

Pi  T 
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The  corresponding  ideal  distribution  of  lift  is  quite  difficult  to  cal¬ 
culate1  and  the  complex  vortex  theory  must  be  employed.  This  has,  in 
fact,  been  done  with  the  aim  of  deriving  a  higher  order  expression  for 
the  lift  distribution  which  could  be  used  in  the  analysis.  Unfortun¬ 
ately  the  final  expression  led  to  extremely  involved  power  integrals 
which  defeated  the  purpose  of  the  study,  therefore  a  simpler  first  order 
approach  was  necessitated.  These  results  are  given  below  without  deri¬ 
vation  in  the  interest  of  brevity: 


Higher  Order  Distribution: 


4C„ 


C,(x)  =  - - - 

n/.25  +  c  /2-cos(tan'1(^/&T/2))\/  x  + 


-  x 


Ct/2 


First  Order  Distributions3  : 


CJt(x) 

C4(*) 


a(x)x  * 

4xCT 

•3a(x) 


for  0.3  s  x  s  1.0 
for  0  s  x  <  0.3 


(16) 


(17) 

(18) 


Figure  5  compares  the  spanwise  lift  distributions  given  by  the  exact 
vortex  theory,  the  higher  order  expression  and  the  first  order  expression, 
respectively.  It  may  be  noted  that  significant  variations  occur  only  in 
the  tip  region.  The  present  analysis  therefore  uses  the  first  order  dis¬ 
tribution  to  calculate  power  which  yields  slightly  optimistic  induced  power 
but  slightly  conservative  profile  power  so  that  the  two  effects  are  essen¬ 
tially  cancelling.  It  is  emphasized  that  the  first  order  distribution  is 
only  a  convenience  for  calculation  and  that  the  correct  distribution  for 
an  actual  minimum  induced  power  design  is  given  by  equation  16.  Present 
conventional  rotors  have  come  close  to  the  minimum  power  assumption  in 

1  simple  momentum  considerations  incorrectly  yield  lift  coefficient  dis¬ 
tributions  proportion  to  l/x  leading  to  finite  tip  loading  and  incorrect 
twist  loading.  This  error  is  most  prominant  at  high  thrust  coefficient. 

2  The  choice  of  the  0.3  station  as  the  point  of  maximum  lift  was  determined 
from  considerations  of  the  maximum  circulation  grading  on  a  wing  of 
finite  span. 
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this  analysis  and  would  come  much  closer  still  if  they  employed  the  non 
linear  twist  implied  hy  equation  16. 

Profile  and  Compressor 

The  concept  of  equivalent  lift-drag  ratio  will  be  used  in  this 

formulation  to  avoid  the  difficulties  involved  with  separate  profile  and 

compressor  power  calculations.  These  difficulties  are  primarily  due  to 

the  effect  of  slot  height  changes.  It  has  been  shown  by  experiment 

(Reference  7)  that  a  change  of  slot  height  for  the  same  momentum  flux, 

or  C  ,  produces  significant  changes  in  lift,  drag  and  compressor  power 
M1 

but  very  little  variation  of  the  parameter  jf/d^  defined  by  Equation  (1). 
An  added  benefit  of  this  method  is  to  provide  a  direct  means  of  relating 
section  efficiency  to  rotor  hover  efficiency  (Figure  of  Merit)  for  any 
airfoil. 

The  equivalent  drag  power  is  defined  as: 


2  3 
PNR  v-3 

p  +  p  =  p  =  - =- 

P  c  \  2 


1 

J 


C/R  C  x3 


dx 


(19) 


in  coefficient  form, 


hi 


C/R  C,x: 

~~WT 


dx 


(20) 


Now  from  the  previous  assumptions  of  uniform  inflow,  the  varia¬ 


tion  of  with  CT  and  x  may  be  written  as: 

C  4ttCT 

*Cl  =  ^T>  0.3  *x*  1.0 


(21) 


c  C  -  x  ^ 
E  ^  ‘  *  .3N 


,  0  s  x  <  .3 


(22) 
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Furthermore,  assuming  a  linear  variation  of  blade  chord  as: 


|  =  *1  -  V  (23) 

where  t±  =  Croqt/R,  t?  =  (CR0()rp  -  C^/R  and  recalling  the  definition 
of  1/ d^  from  Equations  3  and  4  the  profile  and  compressor  power  can  be 
expressed  as  a  function  of  x  and  C  only.  Due  to  the  discontinuity  in 

Ju 

the  definition  of  l/de  the  evaluation  of  the  integral  must  be  accomp¬ 
lished  in  parts  as  shown  schematically  in  Figure  6.  Thus,  the  total 
profile  and  compressor  power  can  be  expressed  as: 


=  C  +  C  +  C  +  C 


(24) 


e  X1  l2  "3 

The  evaluation  of  each  of  these  integrals  is  carried  out  in  Appendix  A. 
Coriolis 


The  last  term  of  the  approximate  analysis  is  the  pumping  or  coriolis 
power  (which  is  required  to  increase  the  angular  momentum  of  the  air  as 
it  flows  out  the  blade  duct  toward  the  tip)  is: 


cor 


=  f2N 


rdF 


cor 


(25) 


where  the  coriolis  force  is  obtained  from  the  equations  of  motion  for  a 
differential  element  of  mass,  dm; 


dF  =  (2uu  x  p)  dm 
cor  K 


or  |dFcor  |  =  2fiVdm 
Where  V  =  p  is  the  duct  velocity. 
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A  change  of  variables  in  made  by  noting: 


Vdm  =  V(r)  p(r)  A(r)  dr  ■  &  (r)  dr, 

where  m(r)  is  the  mass  flux  in  the  blade  so  that  the  mass  flow  rate  of 
the  air  in  the  duct  at  any  radial  station  r  (taking  s  as  an  integration 
variable)  is: 


where  is  the  total  mass  flow  rate  out  of  a  blade 


dm 

dr 


dr 


J 

o 


(26) 


and  —  is  the  mass  flux  per  foot  of  span, 
ar 

The  momentum  coefficient  is  next  defined  generally  by: 


C 

V> 


(r) 


§ 

q(r)c(r) 


(27) 


and  for  hover  by, 


A  small  approximation  is  next  made  in  order  to  solve  the  power  integral 

in  closed  form  by  simply  assuming  V  =  const.  It  should  be  noted  that 

J 

the  assumption  of  constant  V.  implies  that  the  total  head  losses  due  to 

J 

friction  are  identically  balanced  by  the  head  rise  due  to  centrifugal 

concession.  The  power  expended  to  achieve  this  centrifugal  compression 

3 

is  called  coriolis  or  pumping  power. 

A  detailed  internal  flow  analysis  of  a  constant  chord  blade  (Refer¬ 
ence  4)  indicates  that  the  above  assumption  is  approximately  correct  for 
a  ratio  of  Vj/VT  1*2.  In  actuality  the  total  pressure  (hence  jet 
velocity)  first  drops  at  the  inboard  section  and  then  rises  outboard  to 
a  higher  value  than  initially.  Consequently  the  effect  of  assuming  an 

average  value  of  V.  over  the  entire  blade  is  somewhat  pessimistic  but 

J 

certainly  within  the  accuracy  of  the  present  analysis. 

The  mass  efflux  out  the  slot  is: 


da  Pv/R  r 


dx  2V 


c(x)  X  C  (x) 
w* 


And  the  integral  in  Equation  (26)  la  then  as  follows: 


*r  = 


j 

o 


d&  , 

—  dr  = 
dr 


1 

r 


dA  , 
T~  dx 
dx 


J 

o 


*T 


PVT2R 


2V 


j 


1 

/ 


c(x)x  C^(x)dx 


(28) 


3  1 — 

Under  ideal  isentropic  conditions  the  power  expended  for  pumping  is 

regained  by  a  reduced  compressor  power  requirement  (which  is  due  in  turn 
to  the  centrifugal  compression).  The  present  analysis  assumes  that  friction 
losses  prevent  any  increase  in  duct  pressure  and  hence  preclude  any  com¬ 
pressor  power  reduction. 
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By  next  employing  the  original  assumption  that  the  airfoil  is  operating 
at  the  (constant)  angle  of  attack  for  rmim  jj/d,  along  the  blade ,  then 
the  lift  coefficient  can  be  approximated  as: 


C.  =  Ci.C  +  c_ 

i  4  n  5 


(29) 


where:  c^  is  the  lift  curve  slope  with  blowing  coefficient  and  c^  is 
the  lift  contribution  due  to  either  camber  of  angle  of  attack. 


All  of  the  necessary  terms  of  (26)  are  now  available: 


2  2 
pVTTT 


"T  '  2VA 
and  from  (26) 


c  sM 
l  R 


o 


x  dx  -  cr 


J 

o 


(30) 


r 

r 


(fy  - 


PVT2R2 


J 

o 


9/  c(s)C/ 

S2  Ids 


x 

r 


-  Cr 


,2  sill 


ds 


(3D 


The  evaluation  of  the  above  integrals  are  carried  out  in  Appendix  A  with 
the  coriolis  power  coefficient  obtained  as: 


PC0R  pnR2VT3 


Vt 


.0109t1  -  .00262t2 


tl  "  ,3t2 


+0.12399  (32) 


NV  c 

^  (0.itl  -  ,0833t2) 


For  .the  special  case  of  constant  chord,  zero  angle  of  attack  and  zero 
camber  the  above  expression  simplifies  to: 


C 

PC0R 


c4Vj 


(33) 
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Hence  the  pumping  power  can  be  reduced  by  increasing  the  ratio 
c^  =  AC^/AC^,  and  by  increasing  the  ratio  of  Vj/VT  (higher  duct  pressure, 
lower  RFM) .  It  increases  linearly  with  increasing  thrust  coefficient. 

RESULTS 

The  results  will  be  discussed  in  terms  of  the  dimensionless  rotor 
quantities  of  figure  of  merit  and  thrust  coefficient  for  various  values 
of  rotor  taper  ratio  and  solidity.  The  thrust  coefficient-solidity 
ratio,  C^/ct,  will  be  used  to  discuss  the  general  effect  of  blade  loading 
(although  this  parameter  does  not  fully  represent  the  effect  of  varying 
solidity) .  These  dimensionless  quantities  may  easily  be  converted  into 
dimensional  terms  of  power  loading  and  average  blade  loading  by  using 
Figures  4,  7  and  8. 

CORRELATION  WITH  EXPERIMENTAL  DATA 

The  validity  of  the  present  method  can  be  partially  determined  by 
comparison  with  actual  test  data  for  two  NACA  0012  rotors  with  taper 
ratios  of  2  to  1  and  3  to  1  (References  5  and  6).  The  pertinent  charac¬ 
teristics  of  these  rotors  are  listed  below. 


NACA  0012 

NACA  0012 

Taper  ratio 

2:1 

3:1 

Root  chord/radius 

0.253 

0.267 

Tip  chord/radius 

0.127 

0.103 

Twist,  deg 

-12.0 

-8.0 

a 

0.099 

0.095 

Number  of  blades 

2 

2 

The  rotors  were  tested  at  various  tip  Mach  numbers.  As  the  present 
comparison  is  based  on  the  low  speed  section  characteristics  a  correla¬ 
tion  was  attempted  with  the  lowest  tip  speed  case  of  Hpjp  =  0.45.  Only 


one  point  of  comparison  is  provided  by  each  test  corresponding  to  the 
thrust  coefficient  for  which  the  blade  twist  is  given  by  the  moment-urn 
theory  expression  : 


0  (x)  =  Q?i ( x)  -  0^(1. 0) 

where  ar^  is  the  local  induced  angle,  o^(x)  =  tan-1(  Jo^f2/-x) 

Using  the  above  information  the  following  correlations  were  obtained: 


Taper  ratio  2:1  3:1 

Thrust  coefficient  .00378  .00244 

(for  ideal  thrust) 

Experimental  Figure  Merit  0.6l  0.46 

Calculated  Figure  Merit  0.60  0.46 


The  agreement  is  extremely  good  so  that  it  is  felt  that  considerable 
confidence  may  be  placed  in  the  NACA  0012  calculation  as  a  reference 
point  from  which  to  compare  OCR  performance. 

COMPARISON  WITH  DETAILED  ANALYSIS 

The  detailed  numerical  analysis  discussed  previously  was  also  used 
to  determine  the  validity  of  the  present  approach.  This  method  correlates 
well  for  standard  airfoils  and  can  therefore  be  used  with  confidence  for 
the  circulation  control  rotors.  Detailed  CC  section  data  were  input  and 
then  the  rotor  twist  and  slot  distribution  varied  at  each  value  of  C^,  in 
order  to  achieve  the  vortex  theory  ideal  loading.  The  collective  pitch 
and  duct  pressure  were  then  varied  in  order  to  define  the  optimum  figure 
of  merit  at  each  C^,.  The  locus  of  these  points  described  the  limiting 
envelope  of  figure  of  merit  for  the  constant  chord  case  chosen. 

The  results  are  shown  in  Figure  9  where  it  may  be  seen  that,  in 
general,  the  agreement  is  satisfactory  over  the  thrust  coefficient  range 
of  practical  interest  (approximately  .005  -  .015).  For  reference  a 
similar  optimization  has  been  done  which  includes  compressibility  effects 


Again  it  is  noted  that  the  momentum  theory  does  not  yield  the  correct 
twist  distribution  for  minimum  induced  power.  It  was  apparently  used 
for  the  check  case  design  considered  above. 
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obtained  from  Reference  7)  and  is  shown  in  the  figure  also.  No 
attempt  has  been  made  to  calculate  the  compressibility  effect  using  the 
present  analysis  although  it  would  only  require  adjustment  of  the  i/dg 
coefficients. 

EFFECT  OF  SOLIDITY  AND  TAPER  RATIO 

All  results  given  are  for  a  four  bladed  rotor.  The  locus  of 
optimum  figure  of  merit  for  values  of  root  Chord/Radius  =  .03,  .05,  .07, 
and  .09  are  shown  in  Figure  10.  Several  interesting  trends  may  be  noted: 

(i)  The  maximum  figure  of  merit  alwys  occurs  for  the  constant 
chord  case  and  rather  high  0^,/a. 

(ii)  At  lower  C^/o  (which  may  be  the  practical  operating  values) 
the  effect  of  taper  is  beneficial  signifying  important 
related  improvements  in  blade  weight  and  high  speed 
performance . 

(iii)  The  thrust  coefficient  range  for  good  efficiency  is  quite 
large  and  is  maintained  even  at  low  blade  aspect  ratio^. 
This  is  a  significant  result  for  highly  loaded  propellers 
such  as  lift  fans,  ducted  propellers  or  compressors. 

DESIGN  EXAMPLE  AND  COMPARISON  WITH  CONVENTIONAL  ROTOR 

As  a  typical  design  example  consider  a  helicopter  rotor  with  the 
following  features: 

2 

Disc  Loading,  10  pounds/feet 

2 

Blade  Loading,  150  pounds/feet 
Tip  Speed,  600  feet/second 

The  corresponding  dimensionless  values  (from  Figures  7  and  8)  are: 

cT  =  .0117 

CT/e  =  .165 
a  =  .067 

^The  analysis  probably  becomes  increasingly  inaccurate  as  aspect  ratio 
decreases  and  thrust  coefficient  increases  so  that  the  quantitative 
results  under  extreme  conditions  areprobably  optimistic 
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V 


Using  Figure~'10t? ,  it  may  be  noted  that  a  blade  with  CRQqT/R  -  .07 
and  C^p/R  -  .035  will  provide  a  maximum  figure  of  merit  of  0.86  for 
the  CCR.  From  Figure  4  this  is  seen  to  correspond  to  a  power  loading 
of  about  10.2  pounds/horsepower  (S.L.S.),  quite  good  for  a  helicopter 
rotor  system.  Moreover,  it  should  be  noted  that  due  to  the  greater 
section  thickness,  taper  ratio  (2  to  1)  and  more  ideal  load  distribu¬ 
tion  the  rotor  would  be  expected  to  have  a  lower  rotor  weight  fraction 
than  a  conventional  system  . 

A  conventional  rotor  design  would  typically  operate  at  a  lower 

value  of  CT/a  (due  to  blade  stall).  Referring  the  Figure  11  for  the 

NACA  0012  it  may  be  noted  that  two  design  possibilities  arise:  (1) 

For  maximum  efficiency  (FM  -  0.8l)  a  solidity  of  a  =  .0827  is  chosen. 

In  this  case  only  a  small  power  penalty  is  paid  compared  with  the  CCR. 

However,  the  increase  in  solidity  would  probably  result  in  a  consider- 

7 

able  rotor  weight  penalty  ;  (2)  If  solidity  is  reduced  the  weight 
penalty  also  reduces  but  at  the  expense  of  power.  At  a  =  .0756,  the 
figure  of  merit  has  rapidly  dropped  to  FM  =  0.73.  A  comparison  of 
the  two  rotors  for  the  constant  chord  case  (approximately  equal  blade 
weight)  is  shown  for  reference  in  Figure  12. 

Finally  it  should  be  recalled  that  these  results  are  strictly  for 
an  optimum  hover  design.  The  cruise  mode  necessitates  large  reductions 
in  twist  and  variation  in  thickness  ratio  in  order  to  achieve  a  reason¬ 
able  hover-cruise  performance  compromise.  In  general,  for  a  conven¬ 
tional  rotor  this  implies  operation  considerably  off  the  optimum 
described,  herein.  In  contrast  the  CCR  can  still  achieve  close  to  the 
ideal  hover  loading  with  zero  twist  (by  use  of  slot  height  variation) 
so  that  a  significantly  smaller  hover  performance  decrement  would  be 
anticipated,  (Reference  8). 


A  30-percent  thickness  ratio  ellipse  has  10  times  the  flapwise  stiff¬ 
ness  and  7  times  the  torsional  stiffhess  of  a  NACA  0012  section  of 
equal  weight. 

*7 

Rotor  blade  weight  historically  varies  with  the  blade  area  to  the  1.36 
power  while  the  total  rotor  group  weight  historically  varies  in  propor¬ 
tion  to  N1*^( Blade  Wt.)0*95' 
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CONCLUSIONS 


An  integral  solution  for  the  envelope  of  maximum  efficiency  of  a 
Circulation  Control  Rotor  in  hover  has  been  developed  and  shows  good 
agreement  with  a  more  exact  numerical  solution.  The  concept  of  equi¬ 
valent  section  lift-drag  ratio  was  shown  to  be  a  valid  parameter  for 
comparing  blown  and  unblown  airfoils. 

The  Circulation  Control  Rotor  has  very  good  hovering  efficiency 
over  a  large  range  of  rotor  thrust.  The  maximum  values  are  somewhat 
higher  than  the  conventional  NACA  0012  and  occur  over  a  much  wider 
thrust  range. 

Linear  blade  taper  is  beneficial  (compared  to  untapered)  in  the 
useful  rotor  thrust  coefficient  range. 

A  Circulation  Control  Rotor  with  relatively  thick  airfoil  sections, 
linear  taper  and  optimum  spanwise  loading  appears  to  offer  substantia] 
improvements  in  rotor  hover  performance. 


Aviation  and  Surface  Effects  Department 
Naval  Ship  Research  and  Development  Center 
Bethesda,  Maryland  20034 
August  19 ?1 


18 


APPENDIX  A 

EVALUATION  OF  INTEGRALS 


Determination  of  a  closed  form  solution  for  the  profile  and  compres¬ 
sor  power  is  a  lengthy  process  which  includes  the  evaluation  of  four 
integrals.  Equation  A1  is  the  basic  integral  to  be  evaluated  with  the 
limits  of  integration  described  schematically  in  Figure  6. 

(Al) 

where: 


l 


W  cT 

Nx 


for  0.3  £  x  s  1.0 


(A2) 


SCJ =  x 


,  for  0  S  x  <  0.3 


R  =  *1  -  V 

(AU) 

ti  =  croot//r 

(A5'> 

l2  =  ^  CR00T  *  CTIP^/R 

(a6) 

i/d  =  c_C.  for  C  <  C  at  i/d 

(A7) 

max 

*/de  =  ci  +  c?/Ci  for  CA  >  ct  at  */de 

(A8) 

By  substitution  of  Equations  A2  through  A8  in  Equation  Al,  the  power 
integral  can  be  expressed  totally  as  a  function  of  and  x. 
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For  case  I: 


C  <  c’  at  A/d  ;  0  <  x  <  xA 
i  l  e  u 

TQ&X  i  4  J 


where : 


N 


PT  2ttc_ 
1  3 


J 


li  -3 

(-t^x  +  t^x3)  dx 


(A9) 


N 


2trc, 


X 

“1  o. 


t2  5  tl  4 
5"  x  +  IT  x 


(A10) 


For  case  II : 


C.  >  C.  at  Jt/d  ;  x  <  x  *  0.3 
ll  '  e  o, 

max  1 


0.3 


-  .Li2 
p  2tt 

I2 


x 


d  +  ex 


dx 


x  J 
°! 


(All) 


J  = 


UttC^ 

In' 


d  =  c2tx 


4ttC1  , 

6  "  C1  ^ .3N  I  "  C2t2 


KAk 

5d  w 


w  =  d  +  ex 


5vV  ♦  f  v2d2  -  |  »V  *  f 


-,0.3 


d^ln(w) 


(A12) 


For  case  III: 


>  C£  at  l/da  ;  0.3  <  x  i 


max 


x 

r°2 


C 


J 

0.3 


a  +  bx  +  cx 


dx 


(A13) 


where : 


a  =  c. 


b  = 


c^N 


-c2t2N 

c  =  "^cT 


q  =  hac  -o'  q  <  0 


X  =  a  +  bx  +  cx 


V"c#-^ logx  + 


/b2  -  2ac\ 


l 


2c 


.  ^=.tanh 

sh 


_  X 


0.3 


(AlU) 


for  t2  =  0,  constant  chord  blade 
ttC„ 


■I.  b~ 


,5(a  +  bx)  -  2a(a  +  bx)  +  a  log  (a  +  bx) 


(A15) 


0.3 
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For  case  IV: 


C .  <  C .  at  X/d  ;  x  <  x  S  1.0 
t  *  e  ’  o« 

max  d 


N 


p  2ttc, 


N 


1.0 

r 


(-tgX*  +  txx3)dx 


(A16) 


x-' 

°2 


2ttc 


3  L 


-.1.0 


*2  5  \  h 

~5  x  +  T  x 


(A17) 


The  pimping  or  coriolis  power  is  determined  in  a  similar  manner  by- 
evaluating  each  of  the  integrals  separately. 

The  basic  equation  to  be  evaluated  is: 


o 


dm  \ 

ds  J 


ds  = 


pvTV 

2VA 


s2c(s)C. 


r  s  c( s) 


ds  -  c 


J 

o 


5  j  R 
o 


ds  (Al8) 


where  C.  is  defined  by  Equations  (A2)  through  (A6).  Now  the  first  inte- 
ral  on  the  R.H.S.  yields: 


PVT  R 


2VjC4 


s2c(s)CJ 

R 


_  2_2 
pVT  R 

dS  =  ^1  +  I2^ 


(A19) 


where 


4ttC 


(h*  V5  \ 


1  N(.3)2C/R(.3)  \  4  5 


for  0  s  x  <  0.3 


and 


4rrCT  2  (3)^  \ 

I  =  ^-1  i  ir  -  >  'if  -  1  f or  0  s  x  i  1.0 


N  V 


The  second  integral  is  evaluated  as  follows: 


x 

f  2 
s  c(s) 

R 


ds  =  - 


J 

o 


05pVT2R2  ( t.x3 
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Upon  substitution  of  Equation  (A19)  and  (A20)  into  (Al8)  and 
integrating: 
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Cambered  Twenty  Percent  Thiekneee  Circulation 
Control  Ellipse 


NACA  0012  Standard  Reference  Airfoil 


Note:  Chord  Length  for  same  Figure  of  Merit 


Figure  1  -  Airfoil  Section*  used  in  Analysis 
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Figure  2  -  Variation  of  Equivalent  Lift-Drag  Ratio,  for  the 

Circulation  Control  Airfoil 
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Figure  3  -  Variation  of  Equivalent  Lift-Drag  Ratio  for  the 
NACA  0012  Airfoil 
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(Power  Loading)”1,  T/HP 


Figure  b  -  Variation  of  Hover  Power  Loading  with  Disc  Loading 
for  Various  Figures  of  Merit 
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Lift  Distribution 


Non-dimensional  radius,  X 


Figure  5  -  Variation  of  Lift  Distribution  with  Radius  for 
Various  Theories  Blaies) 
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Disk  Loading,  T/S  (lb/ft 
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Figure  8  -  Conversion  Plot  for  Average  Blade  Loading 
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Figure  9  -  Variation  of  Figure  of  Merit  with  Thrust  Coefficient  for 
Various  Theories 


Figure  of  Merit 
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Figure  10  -  Locus  of  Optimum  Figure  of  Merit  for  Circulation  Control 
ROTOR  <«>  W8  ■  -°3 
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Figure  10  (Continued) 

(b)  WR  -  *05 
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Thrust  Coefficient,  CT 
Figure  10  (Continued) 
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Thrust  Coefficient,  CT 
Figure  10  (Concluded) 
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Figure  12  -  Variation  of  Figure  of  Merit  with  Thrust  Coefficient  for  NACA 
0012  and  Circulation  Control  Rotor 
(a)  C/R  -  .03 
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Figure  12  (Continued) 

(b)  C/R  »  .05 
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Figure  of  Merit 


Figure  of  Merit 
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